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Abstract 

We calculate the average persistent current in a mesoscopic metal ring 
threaded by a magnetic flux in the diffusive regime. It is shown that the 
classical electromagnetic energy leads to a universal average current of the 
order of aec/Co, where a is the fine structure constant, —e is the charge of 
the electron, c is the velocity of light, and Co is the classical capacitance of 
the ring. Striking similarities between persistent currents and universal con- 
ductance fluctuations are discovered. We suggest a simple experiment to test 
our theory. 

PACS numbers: 73.50.Bk, 72.10.Bg, 72.15.Rn 
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The experimental confirmation Jl]|| of the existence of persistent currents in mesoscopic 



normal-metal rings threaded by a magnetic flux </> has stimulated many recent theoretical 
works [Q- ||. However, a satisfactory explanation for the experiments has so far not been 
found. In all theoretical models the predicted current is smaller than the experimentally 
measured one. In the present work we shall resolve this discrepancy between theory and 
experiment. We show that in the diffusive regime the energy associated with long- wavelength 
and low-energy charge fluctuations is determined by classical charging energies of suitably 
defined capacitors. The flux dependence of these energies yields the dominant contribution 
to the persistent current. Moreover, we show that the average current is universal precisely 
in the same sense as the average variance of the conductance || . 

This article is divided into two parts. We first study the probability distribution of 
long-wavelength and low-energy density fluctuations in a disordered metal ring in the dif- 
fusive regime. Combining the insight gained from the first part with simple semiclassical 
arguments, we then derive the universal average current and compare it with experiments. 

Consider a thin cylindrical metal ring with circumference L and cross section L\. The 
electrons inside the ring interact with the Coulomb potential, and are elastically scattered 
by impurities. The elastic mean free path t is assumed to satisfy L ^> i > L±. For a proper 
description of Coulomb effects it is crucial to take the neutralizing positive background charge 
density p+(r) of the ions into account. For our purpose it is sufficient to approximate p+(r) 
by a constant p + for r inside the ring, and set p+(r) = if r is outside the ring. Charge 
neutrality requires that the average electronic density equals the density of the positive 
background charge, p(r) = p+(r), where the overbar denotes average over the disorder. In 
terms of the exact retarded/advanced Greens functions G^ A (r, r', e) of spin-cr electrons, the 
true electronic density at temperature T is, for a given realization of the disorder, given by 



We are interested in the probability distribution of the dimensionless random variables 
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p(r)=r oo de/(e)p( 



r, e), where /(e) = [e € ^ T + 1] 1 , and 




(1) 



N, 



k = J dre ike def(e)p(r, e) , k = 1, 2, . . . , (2) 

where the r-integration is over the volume V = LL\ of the ring, and tan# = y/x. We 
have chosen the coordinate system such that the ring lies in the xy-p\ane, with the origin 
at the center of the ring. The cutoff energy E T is large compared with the Thouless energy 
E c = hV/L 2 , but small compared with h/r. Here T> is the diffusion coefficient, and r is the 
elastic lifetime. The N k are the Fourier components of the electronic density due to states 
within an interval of width 2E T around the Fermi energy. Below we show that only these 



states are responsible for the flux dependence of the energy. Because p(r, e) does not depend 
on r, the averages Nk vanish for k ^ 0. More interesting are the variances 



P k = 5N k 5N_ k , (3) 
where SX = X — X for any random variable X. We first calculate 



P(r,r',e,e') = 6p{r,e)5p(r>,e>) . (4) 

The P k can then be obtained from Eqs.|2| and |3]. The essential observation is now that 
P(r, r', e, e') is for |r — r'| ^> i and \e\, \e'\ h/r universal, and can be calculated quite 



reliably within a model of non-interacting electrons ||10|| . The behavior of this function is 
determined by the fact that, at large distances and low energies, density fluctuations in a 
disordered metal are governed by the diffusion equation. As long as inelastic processes do 
not drastically modify the diffusive behavior, their effect can be approximately taken into 
account on a phenomenological level by introducing a finite phase breaking energy in the 
cooperon- and diffuson propagators (see Eq|| below) HlO|fl. To justify this more rigorously, 
let us consider the diagrammatic calculation of P(r, r', e, e'). A general diagram consists 
of two closed loops representing the exact spectral densities p(r, e) and p(r',e'), which are 
connected in all possible ways by static impurity lines. Note that the product of the averages 
is subtracted in Eq.[|, so that the loops must be connected together by at least one impurity 
line. The dominant flux dependent diagram in the absence of interactions is shown in Fig.|l|. 



In the interacting system, Coulomb lines are not permitted to connect the loops. Precisely 
the same diagrammatic rules are encountered in the theory of conductance fluctuations ||. 
The similarity of Fig.|I| with the diagrams shown in Fig. 5 of Lee et al. is evident. Diagrams 
corresponding to interaction corrections to P(r, r', e, e') can be obtained from the diagrams 
shown in Figs. 9 and 10 of Ref. || by simply omitting the current vertices. The contribution 
of these diagrams to P(r, r', e, e') is negligible for exactly the same reasons as discussed in 
Ref. f|. The close connection between fluctuations of the conductance and the density of 
states has also been pointed out in Ref. [|i0| . 

It follows that a reasonable approximation for the function P(r, r', e, e') in the presence 
of interactions can be obtained by performing the calculations for free electrons, and adding 
a phase breaking energy T to the poles of the cooperons and diffusons. The dominant flux 
dependent diagram is shown in Fig.|l[ If e and e' approach the real axis from opposite sides, 
this diagram yields [[Hj 



P(r,r',e,6')=4[^(r,r', e - e ')] 2 , (5) 

where [h 2 / (2nuT 2 )]W(r, r', e) is the usual cooperon propagator, [y is the density of states at 
the Fermi energy for one spin species; the factor of 4 takes into account the spin degeneracy.) 
There exists a similar diagram involving two diffusons, which does not depend on the flux, 
and has been ignored in Eq|| It is understood that the P& calculated below contain only the 
contribution from the cooperon channel. The Fourier transform of W(r, r', e) is for k <C L/t 
and |e| <C E T given by 



i J dr J dr'e ik{9 - e,) W(r, r', e) 



-i 



4n 2 E c (k + 2ip) 2 -ie + T\ , (6) 

where the energy T models inelastic processes, such as Coulomb- and electron-phonon inter- 
actions f§JII|. We have defined the dimensionless flux ip = <f)/<f)o, where 0o = hc/e is the flux 
quantum. The real time Fourier transform of W(r, r', e) is precisely the quasi-probability 
Wt(r, r') introduced in the semiclassical theory of weak localization by Chakravarty and 



Schmid ||11|| . The origin of T is most transparent in this approach: inelastic processes 
destroy the phase coherence between time reversed paths, and are taken into account phe- 
nomenologically by multiplying the quasi-probabilities of the non-interacting system by a 
factor of e~ vt l n . 

Given Eqs.|5] and [], it is straight forward to calculate the coefficients P& defined in Eq.|3]. 
We find 



1 r* " oth|e/(2T)l 



7T 2 m tl^ Jo 4vr 2 P c (m + 2ip) 2 + e + T 
1 
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An 2 E c {m - k + 2ip) 2 + e + T 
Evidently the P& are even periodic functions of ip with period 1/2, so that we may expand 

p(0) oo 

Pk = + E P t ] cos (Anmtp) . (8) 

/ m=l 

After a simple contour integration we obtain 

, 1 r E T /E c j.1/2 



x exp 



-m{x + T/E c ) 1/2 coth [xE c /(2T)] . (9) 



Note that the P fc (0) are for (irk) 2 + T/E c < £ r /£ c proportional to {E T /E C ) 1 / 2 = 0(L/£). The 
cutoff dependence of all non-zero Fourier coefficients P^ m \m > 1, is exponentially small for 
large L/£, so that we may let E T — > oo. In the limit T, F P c , the P^ m ' ) are for m > 1 
universal numbers, determined by the shape of the system. It is important to stress that 
Eq.|9| includes the effect of inelastic processes, which are, however, negligible if the effective 
inelastic energy V is small compared with E c . 

We now develop our semiclassical theory of persistent currents. Consider first the charge 
distribution on the ring associated with the first harmonic N\ of the density fluctuation. 
Although N\ = 0, the typical excess charge on one side of the ring is e(A^ 2 ) 1 / 2 = e(Pi) 1 ^ 2 oc 
e(L/£) 1 / 2 . Because by assumption (L/i) 1 ^ 2 ^> 1, the associated electrostatic energy can be 
calculated classically, and is given by e 2 Pi/(2Ci), where C\ is the capacitance of a thin ring 
consisting of two oppositely charged halfs. The contribution of modes with wavelengths 



large compared with I to the disorder averaged Hartree energy is simply the sum of the 
corresponding classical charging energies, 



A e 2 



^ = E^r^ > ( 10 ) 



k=l 



where the integer A is chosen such that 1 C A < L/f, and the generalized capacitances C^ 



arc 



l l r r ie ik< J>-8') 



r e ik(6-e') 



The sum in Eq.|lO| does not include the term k = 0, because / dr5p(r) = by charge 
neutrality. For k <C L/L± we may approximate to leading logarithmic order Ck ~ Co ~ 
L/[2 ln(L/Lj_)], which is the classical capacitance of the ring. The identification of Eq.pl] 
with classical capacitances is only valid for L±_ <C L. The possibility that charging energies 
could lead to large persistent currents has already been noticed by Imry and Altshuler 
||13|| . Although they examined only the zero mode Co, Imry speculated that local charge 
fluctuations in a globally neutral system might be the key to understand persistent currents. 



Eq.|10] can also be derived formally. A rigorous theorem due to Hohenberg and Kohn [12 
tells us that the exact Hartree energy of all positive and negative charges is a functional of 
the true electronic density p(r), 



e 2 r . f . ,\p(r)-p+] [p{r')-p + ] 



" = 77*7*' ir-V • (12) 

Although Eq.[12| contains explicitly only the bare Coulomb interaction, all higher order in- 
teraction processes are taken into account by demanding that p(r) is, for a given disorder 
potential, the true electronic density of the interacting many-body system. In an infinite 
homogeneous system without disorder Eq.O vanishes, but for a mesoscopic disordered ring 



£1 > 0. Charge neutrality implies p + = p(r), so that the average f2 is uniquely determined 
by the variance of the density. Note the essential role of the positive background charge 
to subtract the product of the average densities. To isolate the flux dependent part of Q, 
it is sufficient to keep only the contribution from modes with wavelengths large compared 
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with £, involving states close to the Fermi energy. According to Eq.[7] only these modes are 
sensitive to the cooperon pole. Using the definitions 0, |3], and [11], we arrive precisely at 
Eq.|T0|. The physical reason for the absence of screening in Eq.[H] is that the random po- 
tential creates in the diffusive regime typical configurations where a large number of excess 
electrons is localized on segments of the ring. Therefore the ring behaves as if it were a 
classical capacitor. 

The sequence of bubble diagrams retained in the random-phase approximation does not 
contribute to the variances Pk- In Fig.||we show the first bubble correction to Fig{l|. Because 
the poles of both Greens functions associated with the loop labeled by e' lie in the same 
half-plane, the r 5 -integration leads to the vanishing of the diagram. Furthermore, if the 
inner left loop in Fig.|2| is not connected by impurity lines to the other loops, the resulting 
diagram is exactly cancelled by another diagram of the same form, where the inner left 
loop represents the density of the positive background charge. This diagram arises from 
the term in the Hamiltonian that represents the interaction between electronic and ionic 
densities. For precisely the same reasons diagrams of this type can be ignored in the theory 
of conductance fluctuations ||. An essential difference between our approach and Ref. |3j 
is that we have included the effect of the positive background from the very beginning in 
a non-perturbative way. The importance of charge neutrality has also been emphasized by 
Schmid [|]. 

From Eqs|||9] and |H] we obtain the Hartree contribution to the average persistent current, 

oo 

7 = -cdU^/dcf) = I (m) sin(47rm^) , (13) 

m=l 

with 

jM = 2afecT & 

7T 2 ^ C k Jo (nmk) 2 + x 2 + 7 m 

xexp [-(x 2 + 7m ) 1/2 ] coth [x 2 /? m /2] , (14) 

where (3 m = E c /(m 2 T), 7 m = m 2 T/E c , and a = e 2 /{Tic) ~ j|= is the fine structure constant. 
The summand in Eq .|14"| vanishes as k~ 2 for large k, so that the value of the sum is dominated 



by the first few terms, and we have set A = oo in Eq.|10[ Eq.[14| is the main result of this 
work. In the limit T, T <C E c , the current is independent of any microscopic details of the 
metal, and depends only on fundamental physical constants and generalized capacitances 
Cki which are completely determined by the large-distance geometry of the system. The 
average variance of the conductance in the diffusive regime is universal in precisely the same 
sense. 

Note that r & E c in the experiments For T ^> E c the current would be exponen- 

tially small. To estimate the magnitude of J*-" 1 -* for 7 m < 1, we may neglect x 2 compared 
with (nmk) 2 in the denominator of Eq.|l4], because the factor e~ x cuts off the integration at 
x = 0(1). In this approximation we obtain 

j(«0 ^ 8 £(2) \n(L/L ± ) c 



9(Pm) , (15) 



I F vr 4 m 2 U7v F 

where I F = (—e)v F /L, ((2) « 1.64, and 

1 r°° r i 

g (/3) = - dxx 2 e~ x coth x 2 p/2 . (16) 

2 Jo L J 

For T <C E c /m 2 we have g{/3 m ) ~ g(oo) = 1. At finite temperature, the relevant energy 
scale where the current begins to deviate from its value at T = is E c . The temperature 
dependence of Eq.|l5| is in agreement with Ambegaokar and Eckern || . As pointed out in Ref . 
|| , this temperature dependence seems to fit the experiment 0. Moreover, 1^ oc vrC 2 , so 
that also the shape of our current agrees with Ref. . However, our current is two orders of 
magnitude larger than the current calculated in all previous works. Assuming L/L± = 100 
and T, T <C E c , we obtain for the first harmonic 1^ ~ —0.9 x I F . The amazing result is 
that has the same order of magnitude as in the ballistic regime M. 

Our theory can be tested experimentally by embedding the ring in an insulator with 
dielectric constant e r . Because capacitances increase linearly with e r , the current should 
then be reduced by a factor of l/e r . The currents measured by Chandrasekhar et al. 
in three different isolated gold rings are of order I F . Although a disorder ensemble is not 
appropriate for a description of this experiment, the fact that the magnitude of the measured 

8 



current agrees with our calculation might indicate that the physics underlying our theory 
is also relevant for the experiment of Ref. @. Levy et al. |]J have measured an average 
current of the order of 0.003 x Ip. The system studied in Ref. [|1|] consists of an array of 
10 7 copper rings with circumference L « 2.2fim, occupying a 7mm 2 area on a sapphire 
substrate. Note that the dielectric constant of sapphire is e r ~ 10, while in the experiment 
of Ref. the substrate consists of Si, which has e r ~ 3. Taking the large dielectric constant 
into account and assuming T = T = 0, our theory predicts for the experiment of Levy et al. 
I m —0.1 x Ip, a factor of 30 larger than experimentally measured. This discrepancy can 
have many explanations. From Eq.|l4| it is clear that an effective r ~ E c can easily explain a 
reduction of / by an order of magnitude. A second possibility is that inter-ring interactions 
and capacitances between rings are important. Note that in the experimental configuration 
of Ref. the average distance between neighboring rings is smaller than L, so that the 
Coulomb interaction between rings is not negligible. 

Our theory is based on the Hartree energy. Quantum mechanical exchange- and corre- 
lation effects have been ignored. We have explicitly verified that the current obtained from 
the first order exchange correction is at least a factor of £/L smaller than the contribution 
of the Hartree term. Because the mechanism leading to the large current has a simple 
semiclassical interpretation, we believe that quantum mechanical correlation effects do not 
essentially modify our result. 

In this paper we have developed a semiclassical theory of persistent currents in the 
diffusive regime. The only quantum mechanical ingredient in our theory are the time reversed 
paths of weak localization. The dominant contribution to the average persistent current is 
universal, and can be derived from the classical electromagnetic energy of the system. 

This work would have never been completed without the advise and criticism of K. B. 
Efetov. I have also profited from discussions with S. Iida, L. P. Levy, and A. Benoit. 
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Figure Captions 



Fig.l: Dominant flux dependent diagram determining P(r, r', e, e'). Solid lines with 
arrow denote disorder averaged Greens functions, dashed lines denote impurity 
scattering, and the maximally crossed ladders define cooperons. 

Fig. 2: First order bubble correction to Fig.Q. The wavy line is the bare Coulomb 
interaction. 
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